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Abstract 

We describe spin transfer in a ferromagnet /normal metal/ferromagnet 
spin- valve point contact. Spin is transferred from the spin-polarized device 
current to the magnetization of the free layer by the mechanism of inco- 
herent magnon emission by electrons. Our approach is based on the rate 
equation for the magnon occupation, using Fermi's golden rule for magnon 
emission and absorption and the non-equilibrium electron distribution for 
a biased spin valve. The magnon emission reduces the magnetization of 
the free layer. For anti-parallel alignment of the magnetizations of the 
layers and at a critical bias a magnon avalanche occurs, characterized by 
a diverging effective magnon temperature. This critical behavior can re- 
sult in magnetization reversal and consequently to suppression of magnon 
emission. However the magnon-magnon scattering can lead to saturation 
of magnon concentration at large but finite value. The further behavior 
depends on the parameters of the system. In particular, gradual evolu- 
tion of magnon concentration followed by the magnetization reversal is 
possible. Another scenario corresponds to step-like increase of magnon 
concentration followed by slow decrease. In the latter case the spike in 
differential resistance is expected due to contribution of electron-magnon 
scattering. A comparison of the obtained results to existing experimental 
data and theoretical approches is given. 
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1 Introduction 



The giant magnetoresistanceT (GMR) of magnetic multilayers is a strong drop 
of the electric resistance of such a multilayer on application of an external mag- 
netic field. This effect arises from the different scattering strength experienced 
by the two spin channels making up the spin-polarized sense current through 
the multilayer and from a change from antiparallel alignment of the magnetiza- 
tion of the ferromagnetic layers to parallel alignment. For antiparallel alignment 
appreciable scattering of electrons occurs in both spin channels, for each spin 
type in those layers where that type is a minority-spin'?. For parallel alignment 
the majority-spin channel experiences negligible scattering, so that it is a low 
resistance channel that short-circuits the high resistance minority-spin channel. 

The dual effect of the GMR is a change of the magnetization state of a 
magnetic layered structure induced by a spin-polarized bias current traversing 
the layers. The mechanism of such a change is transfer of electron spin from 
the current to the magnetization of the layers, as proposed by Slonczewski for 
a ferromagnet/normal metal/ferromagnet (FM/NM/FM) spin valve,3,. This 
is a trilayer structure, with a NM layer sandwiched between two FM layers. 
The spin transfer is predicted^ to lead to steady, coherent precession of the 
magnetization of the layers and, in the presence of a uniaxial anisotropy, to 
switching of the magnetization. These effects should occur for high current 
densities (j — 10^ — 10^ A/cm^) and small lateral dimensions (diameter=100- 
1000 nm). Recently, Katine et o/.0] have reported magnetization switching in 
a spin valve operating in this regime, which they interpret in the framework of 
Slonczewski's theory. However while at small external magnetic fields the be- 
havior was hysteretic just in accordance to Slonszewskii theory at higher mag- 
netic fields the deviations from Slonszewskii predictions were observed. Some 
later, the extended studies of current-driven magnetic evolution were made by 
Urazdin et al. 5 ; the authors reported the hysteretic switching at low external 
magnetic fields and the reversible behavior in high fields. It is also important 
to note that while Slonszewskis theory predict switching-like evolution of mag- 
netization, our studies of a biased nanometer-scale Co/Cu/Co spin-valve point 
contact in] demonstrated a possibility of gradual evolution; similar features were 
observed by Katine et al.0] 

The approach of Ref. is semi-classical: the electron spin is treated quan- 
tum mechanically, but spin transfer is derived from the classical law of angular 
momentum conservation. Further, in Ref. |3| it is ignored that in the spin- 
transfer regime the electron system is strongly out of equilibrium. Actually, for 
the electron system the spin-transfer regime in a spin valve is very similar to 
the regime of point-contact spectroscopy (PCS) of the electron-magnon inter- 
action in ferromagnetic metals. PCS of this interaction was theoretically devel- 
oped by Kulik and Shekhter JOI for homogeneous ferromagnetic point contacts, 
which are short and narrow constrictions between a pair of three-dimensional 
ferromagnetic electrodes. The idea is that in a biased point contact a bias- 
dependent non-equilibrium electron distribution is created. This distribution 
enables energy relaxation of the electrons by incoherent emission of magnons. 
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the elementary excitations of magnetization. Magnon excitation can be detected 
in the second derivative of the current- voltage characteristic of the contact [TT]. 
which is proportional to the magnon density of states. In the same spirit, a 
non-equilibrium electron distribution generated in a current-biased spin valve 
should lead to incoherent emission of magnons in the magnetic layers. To some 
extent Berger ^21 discusses the effect of a non-equilibrium distribution, but his 
intuitive approach lacks a general quantum mechanical basis. 

Further theoretical studies of the non-equilibrium dynamics of magnetic mul- 
tilayers were presented in Refs.^^l where in particular the steps to generalize 
approaches by Slonszewskii and Berger were made. However these studies were 
also based on semiclassical considerations similar to the ones exploited in 

In this article we present a consistent quantum mechanical description of 
spin transfer in a spin-valve point contact, based on Fermi's golden rule and 
taking into account the non-equilibrium electron distribution. Inspired by the 
PCS resultsjjnijd we consider emission and absorption of magnons. How- 
ever, contrary to PCS, which is concerned with the direct effect of electron- 
magnon processes on the electrical transport, we here focus on the effect of 
these processes on the magnetization, which can be strongly reduced by a non- 
equilibrium population of emitted magnons. A change of the magnetization due 
to electron-magnon processes can be probed with the GMR of the contact, as we 
recently demonstrated for a biased nanometer-scale Co/Cu/Co spin- valve point 
contact inj. We will show that different scenarios of the magnetization evolution 
are possible including switching, gradual evolution of magnetization, hysteretic 
and reversible behavior. 

2 Incoherent emission of magnons 

We consider the spin-valve point contact This device comprises two metallic 
electrodes, which are electrically connected through a small circular orifice (di- 
ameter 10-50 nm) in a thin insulating layer. Adjacent to the insulator, em- 
bedded in the left electrode, there is a spin valve with structure FM{ta « 2 
nm)/NM(isp,i ~ ta)/FM{tp « 100 nm), the layer thicknesses being indicated 
in brackets. The thick FM layer is the spin-polarizer. Its fixed saturation mag- 
netization Ms,p points in the positive z-direction. The thin FM layer is the 
spin-analyzer with saturation magnetization q. The direction of M^, ^ is in 
the plane of the layer. It can be parallel or anti-parallel to p. Polarizer 
and analyzer are separated by a NM spacer of thickness tgp.i « 2 nm. This is 
much smaller than the spin-flip diffusion length (and any other scatter length) 
of a NM spacer of a spin valve, so that spin is preserved between polarizer and 
analyzer. The analyzer, in turn, is separated from the insulating layer by a NM 
spacer of thickness tsp,2 {tsp,2 ~ isp,i)- The different scattering experienced by 
the minority- and majority-spin channels in the FM layers is reflected in dif- 
ferent resistivities /9™J{} and p™^/ , which obey Pnm « PpM < Pfm (pnm is 
the NM resistivity). Outside the spin valve the left electrode is continued with 
NM, which is also the material inside the orifice and of the right electrode. The 
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elastic mean free path of the electrons, both in the NM and the FM, is supposed 
small compared to the size of the orifice, so that transport is diffusive on this 
scale. 

The device is biased in the positive direction at a voltage V {V > 0), applying 
—V/2 to the left electrode or reservoir and +V/2 to the right reservoir. The re- 
sulting electron current flowing from the polarizer to the orifice is spin-polarized. 
In zeroth order, i.e. without inelastic processes (inelastic diffusion length ex- 
ceeds orifice diameter) , the resulting spin-dependent electron-distribution func- 
tion in the plane of the orifice is 
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In Eq. (1) the degree of the spin-polarization induced by the polarizer is given by 
a = ARp/Rm = {Pfm ~ Pfm)^pI{t^o^Rm)- This is the difference in polarizer 
resistance seen by minority-spin and majority-spin electrons, normalized to the 
Maxwell resistance Rm , which is the resistance of the corresponding diffusive, 
homogeneous NM point contact, a denotes the electron spin in the polarizer 
(ct = -1-1/2 for majority-spin electrons and a = —1/2 for minority-spin electrons, 
defined with respect to M^^p). Further, /o is the Fcrmi-Dirac distribution, u 
is the total energy of an electron in state k and with spin a, i.e. inclusive the 
electrostatic energy, and e is the elementary charge (e > 0). The distribution 
in the analyzer, since it is so close to the orifice, to a good approximation 
is also given by Eq. (1). The non-equilibrium distribution is similar to that 
of a homogeneous diffusive NM point contact [Tll or mesoscopic diffusive NM 
wires (I5j. It is the average of two Fermi step functions displaced with respect 
to each other by an energy eV, which is the difference of the chemical potentials 
of the reservoirs. In this case, however, the weight of the functions is spin- 
dependent, so that for the energy range of values intermediate between f^ u = 1 
and /k,cr = two values exist: /k.+i/2 = (1 + ck)/2 and /k.-i/2 = 1/2 i for 
majority- and minority-spin electrons, respectively . To arrive at Eq. (1), we 
assume that only the potential drop across the polarizer, because of its large 
thickness and high resistivity, is responsible for the modification of the electron 
distribution with respect to the one corresponding to the homogeneous NM 
point contact (which would correspond to a = in Eq.l). Thus the effect of the 
analyzer is assumed negligible. Further, we concentrate on the spin-dependent 
contribution of the polarizer to the electron distribution of Eq.l, neglecting the 
average over the electron spins. 

Analogous to point-contact spectroscopy of magnons^, ^Jj, the electron 
distribution prepared in a biased spin- valve point contact enables magnon emis- 
sion by electrons in the analyzer, up to a maximum magnon energy eV. Relax- 
ation of created magnons is dominated by absorption by electrons . As usual 
for ferromagnets, we assume that electronic transport is primarily due to sp- 
electrons, so that these electrons control the magnon distribution, irrespective 
the strength of electron-magnon coupling. We further assume that the analyzer 
region exposed to the current is single domain. Finally, at the present level. 
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we neglect escape of created magnons from the region exposed to the current 
and corrections to the distribution given by Eq. (1) due to strong magnon ab- 
sorption by electrons. Thus, applying Fermi's golden rule to magnon emission 
and absorption, and carrying out integrations over initial and final states in 
spherical coordinates, we find the rate equation for the occupation number of 
magnons with energy fkvq in the analyzer: 

^ - ^ / deD{e) J de'D{e') \g\' (1 - (1 + N^)S{e - e' ~ hw^) (2) 

~h',-a (1 - /e,.) N^5{€' ~ e + huoq)]. 

Here D{e) is the electron density of states normalized with respect to the 
unit cell, Ij is an effective matrix element for electron- magnon coupling, i.e. 
renormalized with respect to wave-vector-nonconserving scattering processes 
(see Appendix), and the distribution / is given by Eq. (2). Note that the 
orientation of a in Eq|21 corresponds to minority spins in the analyzer. The 
first term of the integrand applies to emission, the factor (1 -I- N^) denoting the 
sum of spontaneous and stimulated processes. By nature, this magnon emis- 
sion is incoherent, so that incoherent magnetization precessions result. This is 
in contrast to the current-induced coherent magnetization precession described 
in Ref. [3]. Magnons are spin unity quanta. According to spin conservation, 
magnon emission thus is a spin-flip process, which must increase the net spin 
of the electron system by unity by converting analyzer minority-spin electrons 
into majority-spin electrons. Accordingly, the first energy integration in Eq. (2) 
involves minority-spin electrons, while the second integration involves majority- 
spin electrons (e and e' no longer depend on the wave-vector, as angular in- 
tegrations were done already). The second term of the integrand applies to 
absorption, so that the role of minorities and majorities is reversed. When 
]VIs,p and Ms^a are antiparallel, a minority-spin current in the analyzer was a 
majority-spin current in the polarizer, so that for this configuration a = -1-1/2 
in Eq. (1) gives the minority distribution in the analyzer and a = —1/2 the 
majority distribution. For the parallel configuration, minorities and majorities 
conserve their character when they enter the analyzer, so that in Eq. (1) the 
usual sign convention applies. 

To evaluate Eq. (2) for T ~ Q, where the energy of created magnons is 
limited to fiwq < eV , one has to recognize the possible emission and absorp- 
tion processes in the analyzer and to take into account the energy range and 
distribution function (together termed "phase volume") of the electron states 
involved. Magnon emission takes an initial state with fc.+i/2 = (1 + Q^)/2 to a 
final state with _ — 1/2, the phase volume being (1 -t-a)(ey — fiti;q)/4. For 
absorption there is a complementary process from a state with /e'^_i/2 = 1/2 
to a state with = (1 + a)/2, with a phase volume (1 — a){eV — fitJq)/4. 

The sum of these contributions is {a/2)[eV — fiw); the corresponding term give 
rise to a net emission magnons (which corresponds to our choice of the sign of 
a). In addition there are absorption processes from states with f^i -1/2 = 1 to 
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states with /e.+i/2 = (1 + a)/2 and from states with /e',_i/2 = 1/2 to states 
with /e,+i/2 = 0, with phase volumes /ia;q/2 and {l — a)HcOq/2, respectively. The 
total phase volume is fiu}{l — a/2). Note that this contribution to absorption 
exists for purely equilibrium state of the electron system {V = 0, a = 0). In the 
parallel configuration, mutatis mutandis, similar processes with similar phase 
volumes occur. For either configuration Eq. (2) thus goes over into 



dt 



(3) 



= a(Ms_o • M.s,p/2Ms^aMs,p) is a projection of M^^^ on M^^p, weighed by 
the current polarization. Whereas the spin polarization is defined with respect 
to Mg^p, Sz takes into account the sensitivity of magnon-electron processes to 
the polarization with respect to the direction of Mg^o. In Eq. (3) Tm-e ~ 
hr^ \g\ [D{eF)\ twjq is the characteristic time for magnon-electron processes, 
D{eF) being the electron density of states at the Fermi level. 

Eqs. (2) and (3) describe transfer of spin from the spin-polarized current 
to the magnon system. Due to spin transfer the population of non-equilibrium 
magnons increases until a steady state is reached where magnon emission and 
absorption balance each other, i.e where dN^/dt = 0. For magnons of energy 
fujq the steady state is characterized by an effective magnon temperature T^{£, 
obtained by equating the number of such magnons to the average population as 
given by the Planck distribution {ks is Boltzmann's constant): 

reff = 1. e^-^q 

kB 4 1 + (el//fia;q)5/ ^' 

In the limit of weak polarization and for magnon energies fkuq << eV one 
obtains T^^£ « eV/AkB- The effective temperature is larger for Sz < (an- 
tiparallel configuration) than for Sz <0 (parallel configuration). This is natural, 
because the phase volume for magnon creation processes is larger when Sz <0. 
Moreover, in the antiparallel configuration T^l diverges at a critical voltage 
given by 

eVcSz = -fiiOq (5) 

This is interpreted as an unlimited increase of the number of magnons due to 
stimulated emission, resembling a magnon avalanche. This highly excited state 
of the analyzer goes along with a strongly suppressed magnetization and may 
lead to a kind of critical behavior with similarities to the phase transition to 
the normal state at the Curie temperature, which results from strong thermal 
excitation of magnons. As seen in Eq. (4), at the voltages larger than the 
critical one dN^/dt is positive so that has a positive increment. 

Let us introduce the magnon concentration rim normalized with respect to 
elementary cell 
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where fmi^j) is magnon density of states normalized by elementary cell volume. 
It is clear that rim = 1 would correspond to complete suppression of magne- 
tization. Note that a decrease of total magnetization of the ferromagnet |M| 
with an increase of rim for small rim is a well-established fact and describes in 
particular a decrease of |M| with temperature increase. This behavior can be 
interpreted as a result of uncertainty of the orientation of Mg^a within the angle 
^ rim with respect to its orientation at rim = 0. As for the situation of very 
high magnon occupation numbers when n,„ ~ 1, to the best of our knowledge 
it still has not met a consequent theoretical treatment. One could speculate 
that the magnon avalanche leading to rim ~ 1 corresponds to switching since 
for rim — 1 the average, time-independent magnetization of the sample would 
be completely suppressed. At the same time any local fluctuation with the with 
the opposite direction of magnetization would be magnified due to the fact that 
the sign of the projection of the incoming spins on this direction entering the 
driving terms in Eq|2| would correspond to magnon absorption rather than to 
magnon emission. 

However one can expect that the values of rim can be stabilized at some 
rim << 1 due to mechanisms not included in Eq. (4). To look for such mech- 
anisms we should consider magnon kinetics at large occupation numbers (but 
still within the framework of conventional magnon theory implying rim << 1. 

3 Role of magnon-electron processes 

First let us discuss a possible role of magnon-electron interactions. So far we 
assumed that the electron distribution given by Eq.l exists even at high magnon 
concentration. However one can expect that a high magnon emission rate leads 
to decay in the analyzer of the spin-dependent part of the electron distribution. 
The latter consequently can no longer support further strong emission. To 
understand the evolution of this distribution with an increase of the magnon 
concentration we make use of the electron-magnon collision operator: 

/e-™ = j dujv^Vg^'ifeA^ - A', + NA5{e ^ e' - hiu) + NJ{e' -e + hw)] 

-fe',~M - + N^Sie' -£-lTio) + NJ{e - s' + fiiu)]) (6) 

Here i^(ijj) is the magnon density of states. As it is seen, for Nu, >> 1 with a 
neglect of the spontaneous processes this operator tends to establish an electron 
distribution of the sort f minis) = fmajis — fiuj). In addition, the electron- 
magnon processes even at very large N^^ can only take place between of the 
electronic majority and minority states separated by the energy huj since the mi- 
nority electron arising as a result of magnon absorbtion by the majority electron 
can not further absorb magnons and vice versa, the minority electron turning to 
the majority electron by an emission of magnon can not further emit magnons. 
Having these facts in mind one sees that, in particular, the integral of Eq|^is 
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cancelled by the electron distribution 

1 / eV eV \ 

fmaj = 2 [foie + ^ + fi^) + foie - —)j 

fmaj = \ (/o(£ + ^) + /o(e - ^ ~ fi^)) (7) 

So the electron-magnon coupling does not lead to a total decay of spin polariza- 
tion, however the total density of polarized spins supported by the distribution 
resulting from the coupling in question is {al2)fkoD while the initial electron 
distribution given by Eq.[l] corresponds to the density of spins aeV/2D. So one 
concludes that for V > Vc the electron-magnon coupling within the polarizer 
whatever strong it is can not completely suppress the spin polarization. 



4 Role of magnon-magnon processes 

Now let us discuss a possible role of magnon-magnon processes which can be- 
come important at high phonon occupation numbers. Note that these processes 
can be considered as a precursor of the general nonlinearity of the magnon 
physics at high magnon excitation levels. One discriminates between 3- magnon 
processes (originating due dipolar-dipolar interactions and non-conserving total 
spin) and 4-magnon processes related to purely exchange interactions and con- 
serving total spin and total number of magnons (see e.g. |17|'l. For the simplicity 
we will consider purely 3D magnon spectrum {vmil) oc q'^) while generalization 
for the case of more complex spectrum is straigtforward). Let us start from 
3-magnon processes. For the convolution of 2 magnons of the modes 1 and 2 
into a magnon of the mode 3 the efficiency rate can be estimated as ^7] 

l.^(A^,3(l,2)iV,(iV3 + l) (8) 

where Qc is the Curie temperature, M is magnetic moment density while 
73(1, 2) is a dimensionless parameter depending on the details of magnon spec- 
trum. Note that the parameter 7 crucially depends on the details of magnon 
spectrum and on the momentum conservation law. For standard 3D magnon 
spectrum 7 ~ where q2 is the wave vector of the mode 2. However if 
Lo{q — > 0) does not depend on the direction of q momentum conservation law 
allows coalescence processes only if > 3uj{q (see e.g. ^Jj. However in 
general case 

^2 ^ (^^^^ + 2JS{qaf){fiBH, + 2JS{qaf + ATT^ieMsm^ 0^) (9) 

where Hi = Hq — AtthbM + Ha, Oq = ^(QjM) while Hq and Ha are external 
magnetic field and anisotropy field, respectively. As it is seen, in this case 
the restrictions related to momentum coservation are not as severe due to a 
presence of the term depending on 9q. One also notes that while the magnon- 
magnon processes in the bulk of a perfect crystal definitely obey momentum 
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conservation law, it is not necessarily the case for imperfect systems as it is 
shown in particular in the Appendix for electron-magnon interactions. One can 
expect that the efficiency of processes violating the momentum conservation can 
be as high as the for non-violating ones if the spatial scale of inhomogeneity is 
comparable with magnon wavelength. 

In its turn, the rate of 4-magnon process (1,2 3,4) can be estimated as 

- -74(l,2,3)iV2(l + iV3)(l + A^4) (10) 

T4 

where 74(1,2,8) is the corresponding dimensionless parameter which for stan- 
dard 3D spectrum and wave vectors of the same order of magnitude is equal 
to {qa)^ . One notes that since 4-magnon processes simultaneously conserve the 
total number of magnons and the total energy, they can not efficiently modify 
the magnon distribution if it is initially concentrated at lowest possible energies 
(ensuring the largest gain for the magnon avalanche). So we will mainly con- 
centrate on 3-magnon processes. The magnon occupation numbers N can be 
written in terms of our normalized parameters as 

- n^^^ojV-^ (11) 

where V is the phase volume available for the corresponding modes (for standard 
3D spectrum i^(q) V ^ {qa)^)- So one concludes that despite of relatively small 
efficiency of dipolar-dipolar interactions (related to small matrix element /x^M), 
it can be of importance for small magnon wave vectors where the corresponding 
smallness can be compensated even at rim < 1 - by large factors V^^ resulting 
from the phase volume considerations. 

First let us consider a situation when the occupation numbers of the resulting 
magnons N3 are small. In this case one can rewrite Eql^l with an account of 
3-magnon processes as 

—5— = —(S;:eV + fu^)N^-— rim—N^ (12) 

So one sees that the magnon avalanche is stabilized at 

Um^ {\S:,eV\-nw)— — — (13) 

One notes that in this case the gradual evolution of rim is possible with an 
increase of V (starting from Vc) following by the reversal of magnetization at 
some V = Vci corresponding to nm{V) ^1. As for the numerical estimates, 
one notes that for typical ferromagnets Qc/'^'^I^-bM ~ 10'^. Then, according to 
the estimates given in Appendix Tm,e^ ~ 10 — 100. Thus, assuming 73 ~ 10~^ 
one obtains 

n„^(10-3-10-^) ('^-^^'-^) J^ 
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where ij.bM ^ O.OlmV. Thus the values n„i < 1 are still expected for |y — Vc| ~ 
ImV However it may occur that for small values of 73 Vd is too close to Vc to 
allow a wide region for the gradual evolution. 

Now let us consider a situation when the occupation numbers N3 can be 
large which lead to stimulated emission of the corresponding magnons. 

For the mode us one has 

' {S,eV + noj3)N^, 7^-^(1 - 25z) ) + 



dt Tm-ef^3 \ 

27r {^lBM f 73 2 



<(iV.3 + l) (14) 



We will assume that, since lo^ is at least twice larger than uji, SzcV + fko^ > 0. 
So in stationary situation one has 

~ 3^ (15) 



where 



while 



So if /? < the stationary solution for N^, corresponds to 

{\S,eV\-hij)-——-^ (l_I^)=n„ (18) 



2 2 / ^ -2 \ 1/2 



or 

^2 r|2 / 

As it is seen, the behavior depends on the ratio 

hm ~ (|5,ey-MI ey^(r„_ea;)i/2^^ ^ ' 

Namely, in the first scenario hm < ric 

rim ~ fim (21) 

In this case the gradual evolution of Um is possible provided Vd — Vc is large 
enough. In the second scenario hm > tie one has 

Um ~ ric (22) 
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While at V — Vc ^ one always has n„i ~ n„i, at higher V the condition 
fim > n-c can be reached at some V — Vc2 provided Vc2 < Vd- It is important 
to note that when rim ~ n-c the value of rim decreases with an increase of 
since the magnon generation corresponds to SzeV < 0. 

Thus, the non-monotonous evolution of magnetization is expected: at critical 
bias SzcV ~ —hui the value of rim — = Vc) is reached while the further 

increase of V leads to a decrease of the magnons number. One notes that in 
this regime 



and thus N^^ increases with the increase of V . The situation again becomes 
critical when the values of becomes comparable with N and a situation of 
energy diffusion is established. However we will not discuss this situation in 
more detail. 

Note that the actual physical picture with an account of magnon-magnon 
processes is rather complex and sensitive to details of magnon spectrum etc. 
Thus our simplified analysis aims only to reveal the main features. We can 
conclude that the magnon-magnon processes can lead to saturation of magnon 
emissions at rim < I7 that when no switching occurs. In this case the fur- 
ther gradual evolution of magnetization with an increase of V is possible. The 
scenario depends on the relation between ric and fim- 

5 Comparison to previous theoretical approaches 

A short comparison with the approach of Slonczewski^j is in place. In the latter 
case the spin transfer from the incident electrons to the layer is proportional 
to the amplitude of magnetization precession and is zero for M being parallel 
or antiparallel to Mq while in our case the spontaneous processes exists as well 
allowing spin relaxation even for the spins parallel or antiparallel to Mq. In the 
case of mechanism by Slonczewskii the spin transfer from the incident electrons is 
equal to js,inj sinfl provided d> Ip ^ {h/pp)eF I E^x where Eg^ is the exchange 
energy, Ip is the spin precession length and js,inj is the spin current. Since E^x is 
not too small with respect to the Fermi energy this condition holds for practical 
values of ta- Thus this mechanism is restricted by the near-surface layer with 
a thickness Ip. The evolution of the precession angle is described by Eq.l7, the 
key equation of the paper, which can be in our terms be written as 



where a is the Gilbert parameter while js,inj is the spin current in the analyzer. 
(Note, that, strictly speaking, the apphcability of Gilbert damping, suggesting 
the equilibrium state of the ferromagnet, is questionable for the non-equilibtium 
electron distribution). To compare with our considerations we first note that 
for small 9 the number of " coherently excited magnons" can be related to 9 as 
^m ^ 9"^ 12. Then, one can consider the product auj as the magnon relaxation 



d0_ 
'At 




)sa\9 



(23) 
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rate, l/rm- We believe that the dominant mechanism of the magnon relaxation 
is related to magnon-electron coupling discussed above. Thus one can rewrite 
this equation as 

^ = - (24) 

As it is seen, the criterion for the current-driven evolution of magnetization can 
be written for the Slonczewskii mechanism as 

js,inj > 7 (25) 
Tm) 

which can be rewritten in terms of the bias with an account of the estimate of 
Tm given above as 

V>r]-'^Vc (26) 

where rj = js,inj /jg'^inj while j^""^ is the spin current which would exist under 
the same bias provided the structure would be ballistic. One notes that this 
criterion differs from that given by EqOby a factor > 1 in the r.h.s. of the 
inequality which for the diffusive electron transport is much larger than unity. 
Thus there is a broad region of biases 

huj < \Sz\eV < rj^^huj 

where the excitation of magnons according to our scenario is possible while the 
mechanism by |3] is not efficient. 

Note however that the approach suggested in our paper is completely self- 
consistent one and is limited by quantum-mechanical treatment while the ap- 
proach of [3] exploits a combination of quantum-mechanical treatment of elec- 
tron spin evolution and semi-classical transfer of x— component of electron spin 
to the layer as a whole described with a help of classical Landau-Lifshits equa- 
tion. As a result, some questions do not meet unequivocal treatment - in partic- 
ular a mechanism leading to a change of z— component of spin by portions equal 
(1/2) which governs the spin evolution in a consequent quantum-mechanical pic- 
ture and is expected to have "non-coherent" character. It is also not clear in 
what way the single-electron spin is transferred to the excited region as a whole 
since any given electron is coupled only to its closest surrounding. 

To some extent these questions were noted by Berger ^21 who exploited the 
similar mechanism of boundary-induced spin transfer to describe emission of 
spin waves by incoming electron spin flux. He attempted to write the equation 
for z— component of electron spin, but he simply postulated it in the most simple 
relaxation time approximation form to fit a quantum-mechanical principle that 
z-component of spin can be changed by portions of 1/2 (Eqs. (9-12) in 52]) ■ In 
this way he also came to a conclusion concerning instability of magnon system for 
the current exceeding some threshold value similar to one we have formulated 
on the base of our consistent quantum- mechanical treatment. Note however 
that, as the paper this approach does not take into account a possibility of 
non-equilibrium electron distribution, relating the chemical potential difference 
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majority and minority subbands to a difference in pumping of electron momenta 
at the boundary due to a difference in drift currents. To our opinion such an 
approach is equivalent to an assumption that the energy diffusion length is 
equal to elastic mean free path and strongly underestimates the real potential 
difference. 

Since one notes similarities between our results and results of jS] , includ- 
ing a possibility of switching discussed above, magnon escape factor etc. one 
can consider our scheme as a generalization of approaches by Slonzcevskii and 
Berger. But we also predict a possibility of stabilization in magnon system at 
the levels still much less than corresponding to switching - in contrast to results 
by 0, [El. 

As for high excitation levels V — Vc > 1 we would like to note that there is 
still a problem whether the real magnetization evolution is in accordance to our 
"incoherent" scenario or it has features of "coherent" classical magnetization 
evolution of the sort considered in (with some correction of the Gilbert pa- 
rameter a for the nonequilibrium case) . In any case one expects that if initially 
M and Mq are parallel or antiparallel an initial stage is described by our scenario 
since it is related to magnon excitation which can be considered as quantum 
fluctuations and exist even for 6 = 0. The further behavior is not as clear since, 
as we have noted above, the situation of very high magnon occupation numbers 
is still not perfectly clear. In particular one can expect a sort of "producing 
of coherency" originated due to effective magnon-magnon processes and estab- 
lishing the "coherent" evolution. Note that a similar problem was discussed 
for acoustoelectric generators of sound waves where the coherent signal arose 
(or not arose) from initial electron-drift driven emission of incoherent phonons 
(CHI)- Certainly this choice can crucially dependent on the factor of boundaries 
which can support - or destroy - such coherency, and on the excitation level 
(that is on the product \Sz\eV). 

We would like also to emphasize that for non-ballistic limit where the electron- 
magnon coupling in our case has a bulk character (see Appendix) in contrast 
to mechanisms of spin transfer considered in [3], ^21 our "bulk" mechanism is 
expected to dominate in any case. 

To conclude this Section, we would like in addition to consider a situation 
when the exciting current passes only through small area of the analyzing layer 
which takes place when the planar multilayered structure is excited with a help 
of a point contact attached to such a structure. This geometry corresponds in 
particular to experimental situation of the papers [J], [H] which will be to some 
extent discussed in what follows. In contrast to the picture suggested above, 
the emitted magnons have a possibility to escape the excited region propagating 
along the plane of the analyzing layer (assumed to have an infinite area). If one 
demotes a size of the excited region (being equal to the size of the point contact) 
as a, the escape time t^"^ is expected to be equal to a/vm where Vm — hkjn/rrijn', 
here Vm is a magnon velocity while fc,„ is the magnon wave vector with a smallest 
possible value of ~ 1/a. Correspondingly, the largest possible value for t^'^ is 
obviously a'^nim/h. 

This escape can be taken into account by adding to Eq. (3) the relax- 
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ation term — TV^^/r^f'^, where r^f'^ is the magnon escape time for the point- 
contact geometry. This gives Tm-e/T^'^ as extra term in the denominator of 
Eq. (4). Divergence of T^/£ is not prevented by magnon escape, but fast escape 
{'Tm-e/T^'^ >> 1) Icads to 'A strong increase of the critical voltage according to 

eV^S, « ^(7:,„_,/r^^^)fia;q (27) 

Taking into account that r^ig ~ ujUkpd) one obtains 

r™-e/rr > (28) 

where the equality corresponds to the smallest magnon in-plane wave vector 
~ 1/a. As it is seen, in this case the criticality starts at the lowest possible 
in-plane wave vector. So the criticality criterion is in principle independent of 
the magnon frequency provided the ratio given by Eg 1281 is larger than unity. 

Note that according to Eq[57|in the corresponding regime Vc cc l/(a^). For 
diffusive electron transport the critical current is Ic = Vc/Rm is thus inde- 
pendent of the orifice size which agrees with prediction by Slonczewski for the 
similar geometry although in terms of coherent excitation of spin waves 



6 Comparison to experimental situation 

Thus we have shown that the incoherent stimulated emission of magnons 
trigger an avalanche evolution of magnetization at biases smaller than ones 
necessary for the switching behavior according to the mechanism suggested by 
Slonczewski. However the final stage at these moderate biases corresponds to the 
state with very high magnon temperature but still with the initial direction of 
magnetization. Such a behavior is different from "switching" predicted in In 
our scenario the first step-like evolution of magnetization with the applied bias 
(at V = Vc = \hujjn/eSz\ can be followed by the gradual increase of the number 
of magnons and, correspondingly, by a gradual decrease of the z-component of 
total magnetization. Note that neither the approach by Slonczewski nor the 
approach by Berger have predicted such a gradual evolution. Unfortunately 
the actual physical picture at this stage which is controlled by magnon-magnon 
processes is rather complex and sensitive to details of magnon spectrum etc. 
Thus our simplified analysis allows only to reveal the main features and does 
not allow direct quantitative comparison to experimental data. 

First we note that qualitatively many experiments on current-driven excita- 
tions in magnetic multilayers |ZI,[H|i ^ilHI; lEI unambiguously exhibited signa- 
tures of a gradual evolution of magnetic state with an increase of bias. 

In the papers [Sj, ^ current-induced switching of magnetic moments in 
Co/Cu/Co pillar structures was reported. The switching exhibited hysteretic 
behavior as function of injected current and external magnetic field and was 
interpreted in terms of theory by Slonczewski ^2. • Note however that the initial 
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phase of the process clearly indicated a gradual evolution with an increase of the 
current following step-like increase of differential resistance. At high external 
magnetic fields instead of switching the data exhibited spikes for one of the 
directions of the current which were attributed to emission of spin waves. Since 
the theory [S] exploiting a constant value of Gilbert parameter does not allow 
an existence of stable spin-precessing mode, the authors of 4 explained the 
observed behavior as a result of a dependence of a on the angle of precession. 
Although the latter mechanism can in principle explain spike-like features, it 
can hardly explains the gradual evolution. Note that the gradual evolution 
of contact resistance with an increase of the current depending on the current 
direction was also observed in 

The transition from hysteretic behavior to the spike- like one with an increase 
of magnetic field was also reported in |5] for permalloy-based nanopillars. We 
would like to note that authors of have suggested a possible role of nonequi- 
librium magnons with high temperatures to explain the behavior observed. The 
also observed a presence of telegraph noise in the vicinity of the spike. It is 
important to note that the noise was pronounced only at very narrow region of 
biases and has not revealed any trace of hysteretic behavior. This fact could 
not be easily understood suggesting magnetization reversal since in the latter 
case the anisotropy field would give rise to a hysteretic behavior. 

To our opinion, the spike can be related to stabilization of the magnon 
occupation numbers due to magnon-magnon scattering according to the scenario 
(2) considered above. Indeed, for the resulting "incoherent precession state" 
the direction of magnetization is only slightly "smeared" with respect to its 
original direction and thus no hysteretic behavior is expected. In contrast, 
the variation of the bias around the corresponding threshold value leads to a 
reversible behavior of the magnon occupation number. The step-like increase 
of resistance leading to the spike in differential resistance can be related to 
additional resistance related to effective electron-magnon scattering rather than 
to magnetization reversal. 

Correspondingly, the transition from hysteretic behavior to spike-like behav- 
ior can be related to the transition from scenario (1) to scenario (2). Indeed, 
with an increase of magnetic field a role of the angle-dependent term in [^de- 
creases with respect to the term ^sHi. As a result, the the phase volume for 
3-magnon processes decreases thus leading to a decrease of 73 . As it is seen, it 
is favorable for regime (2) where an increase of bias does not lead to an increase 
of Um at least in some interval of biases. 

In the papers [7], [HI magnetic mulitilayers wee excited with a help of point 
contact. Correspondingly, the excitation took place only at small area of the 
multilayer concentrated around the point contact which is in agreement with 
model discussed at the end of our Sec. 2. The differential resistance spectra 
obtained in [7| for Co/Cu multilayers excited by the metallic point contact 
clearly demonstrated a gradual increase of resistance following an initial spike 
which the authors ascribed to stimulated emission of magnons according to 
scenario [T^ . 

To our opinion, the experimental data mentioned above evidence the inco- 
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herent character of magnons emission (at least at the initial stage) since it is 
the latter mechanism that can - according to our model - explain co-existence 
of "critical" spike-like and" switching" behavior with a gradual evolution. 



7 Conclusions 

To conclude, we have given a consequent quantum-mechanical treatment of 
spin-polarized transport in spin-valve point contact with an account of non- 
equilibrium electron distribution is given. It is shown that at large biases 
an avalanchc-likc creation of low-frequency magnons within the ferromagnetic 
layer up to is possible which agrees with earlier predictions based on semi- 
phenomenological models. However in contrast to earlier studies [5] it is found 
that the stabilization within the magnon system can be achieved at finite (al- 
though large) magnon temperatures and the gradual evolution of such a state 
with the bias is possible. These results are in agreement with existing experi- 
mental data. 



8 Appendix 

Let us consider an effect of boundaries and disorder on the magnon-electron 
matrix elements. Following Mills et al. ^Hl we will write the element corre- 
sponding to a transition from the electron state k to the electron state k' in a 
presence of created magnon in the state q as 

A1|k^k = ./^^ y dr3expi(k-k'-q)r (29) 

where s and M are the electron spin and sample magnetization, J is an exchange 
constant, a is the lattice constant, k, k' and q are the wave vectors of initial and 
final states of electrons and the magnon wave vector, respectively, while V is the 
normalizing volume. One notes that typically the integration gives a standard 
momentum conservation law. However now we are interested in a case when 
the momentum conservation is violated which in particular is related to a finite 
size of the layer with a thickness ta in z-direction. Thus we will assume that in 
plane of the layer the momentum conservation holds in a standard way and will 
concentrate on the integration over x. The corresponding factor arising in the 
expression for |A^|^ can be readily written as 

(30) 



For a given initial and final electron energies s — £k,-(T and s' — Ek'.o- One 
notes that the Fermi surfaces are typically separated by relatively large gap 
Akp — kpEex/sp- For long- wave magnons we expect q << A.kp which allows 
us to neglect q in our estimates. We will first integrate over angular variables 
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t?, if, Lp' . Having in mind that the difference e — e' (controlled by the distribu- 
tion functions) is much less than E^x we also will neglect this difference in course 
of the angular integration. Thus one has kx = kp- cos ■(9 and k'^ = kF,+ cost?'. 
The momentum conservation in perpendicular direction eliminates integrations 
over (f' and "d' with an obvious result sint?' = {kF-/kF,+) sai9 (which is related 
to a conservation of the in-plane component of the momentum). Correspond- 
ingly, one obtains 



cosi9' = (1 - {kF,-/kF,+fsin^ 4) 



1/2 



One clearly see a gap preventing small values of cos??' (since the modulus of 
x-component of electron wave vector for majority electron is larger than for 
minority) and finally obtains for the denominator in Ea l3UI 

V kF e 



tklU-Jeil-2^)^2^] ^tlkli^i^-^] (31) 



where ^ = cost?. Thus performing the integration over ^ one has a factor 

a 



tlAkl 



Integrating over angular variables corresponding to the initial and final states 
one notes that the effective matrix element can be estimated as 

|7W|' = J^{kF/{AkFfta) (32) 

The effective rate equation can be rewritten as 

^ = ^E /de^e) y'd£'K£')l-MP(/e..(l-A',-.)(l+AC)<5(£-e' + M- 

-(1 - fe,a)fe',-M{e - s' - hw^)) = (33) 

Here v is electron density of states normalized with respect to elementary cell. 
One notes that in a view of the normalizing factors of electron wave functions we 
have used in course of matrix clement calculation one should also put ^ ^a, 
X^k' ^iid corresponding factors eliminate the normalizing factor 

arising in course of our estimates of the matrix element which is taken into 
account in Eg 1331 

Note than if one would write the similar equation (integrated over the angu- 
lar variables) in the case when momentum conservation over x is allowed (when 
AkF = the resulting value of denominator in Eg 1301 would be limited only by 
the smallest possible wave vector along x that is kx — T^/ta and as a result of 
integration over ^ one would have 

a 
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Correspondingly, an integration over angular variables would give a factor 

kp 

q 

because in the limit g — > in the elastic limit one has a divergence of the 
corresponding cross-section. Thus the matrix element within the equation 1331 
would have a form 

\M\^ - (34) 

that is a decrease of the magnon-electron rate with respect to a case when 
momentum conservation is obeyed is descrimed by a factor 



One notes that the finite electron-magnon coupling strength for low magnon 
frequencies (and relatively small biases) is in our case achieved due to boundary- 
induced momentum non-conservation. One expects qualitatively similar effect 
from any factors breaking momentum conservation. The effect of alloy disor- 
der significantly increasing an efhciency of magnon-electron processes in low- 
frequency region was considered in although in this case the disorder was 
related to on-site exchange constant. We believe that the elastic scattering of 
electrons can also enhance the coupling. 

Let us first consider a role of this scattering when the analyzer is still thin 
with respect to the elastic mean free path le- One notes that a presence of 
elastic scatterers leads to a modification of electron wave functions with respect 
to plane waves, due to the scattering. Each of the scatterers produces the 
scattered wave, so the plane wave in the lowest approximation with respect to 
the scattering processes (which holds when ta < le) as 

»kr »kr , f expzkRi/c|Rir 

e — > e t x + 



|Rir 

where i numerates the scatterers, Ri is the coordinate of a scatterer while / is 
the scattering amplitude. One notes that the replacement of the plane waves in 
Ea l29l bv the modified ones gives a non-zero result even with no account of the 
boundaries. Since the magnon wave vector is assumed to be much less than the 
electron ones the integrand in the Eql^Hlfor each of the scatterer has a form 

gikr , g-i(kRi + fe|rRi|) ^ gi(kr+fcr) , ^iAkr , gi(kk)Ri ^gg-j 

where r = rRi. With an assumption Afc << A: the integration over r gives 

expi(kk)Ri (37) 



jfcAfc 



where \A\ ^ 1. The total matrix element is a sum over all the scatterers in the 
sample. However in the square of this element the products of contributions 
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of different scatterers vanish due to strongly oscillating factor in Eq|^ depend- 
ing on the scatterer coordinate. As a result, one finally obtains the following 
estimate: 

\M\'-^ (38) 

For this estimate we have taken into account that f'^Ni ^ and that ~ 
e^ikp/^kf- As it is seen, for ta < le the contribution of scattering is smaller 
than the effect of boundaries. However one expects that the result of Eal38l 
holds even if for the diffusive transport in the analyzer. Indeed, one notes that 
actually the contribution of each of the scatterer is formed at distances ^ l/(Afc) 
from the scatterer. So the derivation implies a ballistic transport only at this 
spatial scale which is assumed to be less than ta- 

Note that in this case the electron-induced magnon relaxation rate with an 
account of the fact that 5k ^ kpiJ /sp) can be estimated as 

^ kk^ 

and thus the factor (l/kple) plays a role of Gilbert damping parameter. 
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